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Abstract. We calculate the boundary correlation function of fixed-to-free boundary 
condition changing operators in the square-lattice Ising model. The correlation 
function is expressed in four different ways using 2x2 block Toeplitz determinants. We 
show that these can be transformed into a scalar Toeplitz determinant when the size 
of the matrix is even. To know the asymptotic behavior of the correlation function at 
large distance we calculate the asymptotic behavior of this scalar Toeplitz determinant 
using the Szego's theorem and the Fisher-Hartwig theorem. At the critical temperature 
we confirm the power-law behavior of the correlation function predicted by conformal 
field theory. 
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1. History and outline 

Onsager's solution pQ for the two dimensional Ising model proved to be an important 
stepping stone for the development of a theory of critical phenomena, and for the 
invention of conformal field theory (CFT). The solution has been studied by various 
methods: by using essentially the transfer matrix technique and the Jordan- Wigner 
transformation [21 [3j HJ E], by taking a combinatorial approach of using the Dimer 
problem [6j U\ E], by using Grassman variables [9], and so on. These methods lead 
to many exact results including the evaluation of the spin-spin correlation function 
[HIS, [TO], the boundary magnetization and the boundary spin-spin correlation functions 
[8], the interface free energy and the magnetization profile across the interface [5], and 
so on. 

At the critical temperature many of these results are predicted by conformal field 
theory (CFT). The rational CFT [TT] of the central charge c = | describes the Ising 
system at the critical temperature, producing correct operator contents and correlation 
functions [12]. It also describes the conformally invariant boundary conditions and 
predicts the dimensions of boundary operators [13] • 

In this paper we will calculate the correlation function of fixed-to-free boundary- 
condition-changing (bcc) operator directly from the lattice model. By doing so we 
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will confirm one of the predictions of CFT, namely, the dimension of fixed-to-free bcc 
operator, and also obtain the off-critical behavior of the correlation function. We use the 
formulation by Kadanoff pi] to write the general boundary correlation functions £ls ct 2 x 2 
block Toeplitz determinant. The asymptotic behavior of the Toeplitz determinant when 
the size of the matrix goes to infinity is, in many cases, given by the Szego's theorem 
[H]. At the critical temperature, however, the Szego's theorem does not capture the 
asymptotic behavior of the determinants and we must use the Fisher-Hartwig theorem 
to deal with the singularity. This is possible since we can transform our block Toeplitz 
determinant into a scalar Toeplitz determinant when the size of the matrix is even. 
The basic definitions and theorems about Toeplitz system are collected in Appendix. 



2. Boundary correlation function: General 

2.1. Partition function 

We consider a 2-dimensional Ising model on a square lattice. Lattice sites are labeled 
by two integers (j, k) representing x and y coordinates. The spin variable at (j, k) is 
denoted by Sjk = ±1. The partition function of the Ising system is given by 



j,k=—oo 



(1) 



z = ex p 

W 

where the horizontal coupling, Kj^, is assigned to the bond connecting (j, k) and k) 
and the vertical coupling, Kj^, is assigned to the bond connecting (j, k) and (j, k + 1). 
For the horizontal coupling, Kj t k, we also define the dual coupling, 9j t k, such that 
tanh^fc = exp(— 2Kj tk ). 

It is derived in [I] that the above quantity is evaluated as 



IJ 2 sinh 2K jtk y/det [Q - e"^] , (2) 
vJ'.* / 
where Q and e~ iPx are infinite matrices that we describe next. 

These matrices have a block structure: the components of the matrix are labeled by 
indices j, k and r where the last one takes either 1 or 2j|| First we define the translation 
matrices as 

[e lPx )j,k,T,j',k',T' — $k,k'0~j±l,j'0~T,T', 

[e 1PV ]j,k,T,j',k',T' — 0~k±l,k>0~j,j'0~T,T'- 

These simply raise or lower the corresponding indices by one. We also define partial- 
translation matrices, T±, such that they raise or lower the composite index (2k + r) by 
one. Using 2x2 notation to explicitly show their action on r-space, they are written 
as T± = cxi exp(±iY) ex p(i^-o"3) or, equivalently, as 

V i o J + v eiPv / 

| This reflects the two types of fermions in the model [I]. 
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using the standard Pauli matrices, 

* 1= (i I) a2= {°i ~o) a3= {l -1 )• 

Using these (partial-)translation matrices we now define the Q-matrix as 

Q = T + e~ 2@a2 T_ e^e-^^e"^ , (3) 
where and K are diagonal matrices defined by the (dual) couplings as 

Qj,k,T,j',k',T> = 0~j,j>5k,k>5 T ,T' Qk,j 
Kj,fc,Tj',fc',r' = &j,j'&k,k'&T,T'Kk,j ■ 

The partition function ([2]) is invariant under the following similarity transformation. 

Q _> Q = e - 2@a2 T_ e ip *e- 2K<T2 e- ip * T + . (4) 

2.2. Boundary correlation function 

Using (j2J) we will derive the expectation value of two boundary operators separated by 
n bond-lengths. First, following [4], let us briefly describe how the expectation value of 
an operator in the bulk may be evaluated. 

An operator O may be realized on the lattice by changing the couplings in a specific 
way (see, for instance, [15J). In such case the insertion of O will modify Q, say, to Q'. 
Then, (O) may be written using (j2J) as 



Z[Q] \\ \ W s inh2i^ jfe 1 ' ( 



Q'-Q 

Q - e-w* 



(5) 



where we use the notation Z[Q] to emphasize the coupling dependence of Z. If Q' and 
Q differ only by a finite number of couplings then the size of the matrix that one needs 
to take the determinant is finite and given by the rank of Q' — Q. 

Let us define Qo as a Q-matrix for the uniform but possibly non-isotropic system, 
i.e., (Kj t k, Kj,k) = (K, K) for all j and k. 

Qo = exp(-if a 3 ) e 29 ^ exp(if a 3 ) e~ 2K °\ (6) 

where 6 is the dual coupling of K. We can decompose Qo as 

Q = An + A- 1 (i-n), (7) 

where the eigenvalue A and the projection II are written as 



TT - - 4. P K ^rr c cosp y -c-ia 3 s sin p Kn 

2 + 2 2v^T ' [) 

using the following convenient abbreviations: 

(c, s,c,s) = (cosh 2K, sinh 2K, coth 2K , csch 2K) , 
£ = cc — ss cospy . 
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Note that the bar means the horizontal bonds, not complex conjugation. One can do 
the same for Q (J3J) and define Qo and II in a natural way. A remains the same. 

The decomposition (I7j) is useful when one evaluates, for instance, a matrix inversion 

(Qo- 



-ipaA— 1 



Qc 




n 



+ 



n 



for m < 
for m > 



For m = the above becomes 
"(Qo-e-^)" 1 



A _1 n, 



jj 



(9) 



which will be useful shortly. 

Similarly to adding an operator, one can also add a boundary to the system by 
changing the couplings in a specific way. Let us add a boundary vertically at j = 0, 
producing two half-planes on both sides of the boundary. For examples, setting K ^ = 
for all k will impose a free boundary condition, and setting Kq^ = oo will impose a 
fixed boundary condition. We define Qb as the corresponding Q-matrix. By definition 
the difference Qb — Qo has non-zero components only at the diagonal block of j = 0. 

To evaluate (O)b, the expectation value of O in the presence of the boundary, we 
define Q' by implementing the insertion of O into Qb- Then, (O)b is written as the 
following ratio 



sinh 2K 



in) ML 1 1 

[ )h ~Z[Q b ] \iisinh2^ )fe 



det 



l + (Qc 



>-ipx\-l 



(Q' - Qo) 



i + (Q c 



HQb-Qo) 



(10) 



where K h stands for the horizontal couplings as appeared in Q b , and K 1 stands for those 
in Q'. 

To simplify the above equation let us define a Laurent matrix (see Appendix for 
the definition) Ab and a projection r\ x as to satisfy 



Qb - Qo = A b Vx and h 



X}j,k,T,j',k',T' 







for j = 
otherwise 



The determinant in (1101) may be rewritten as 

1 



det 



1 + 



1 



1 + (Qo - e-^)- 1 A h rj x Q - e"^ 
Using (jSJ) we obtain the identity 

1 



l + (Qo-e- J f-)- 1 A b % 
and the determinant ( flTT) becomes 

1 



l-(Qc 



(Q' - Qb) 



1 A b rj x 



det 



1+1 



Qo - e iPx 



A b ?7 2 



1 



l + A-inAb/ Qo-e-^ 



l + A-^Ab 
(Q 7 - Qb) 
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Now defining a scalar Ab as to satisfy 

(A 6 -A)n = nA b n, 

the determinant in ( fTUj) finally becomes 



det 



1 + 



Qo - e-to» 



l-A bVx (l-Ab'lIAb 



Qo 



(Q' - Q b ) 



(12) 



Without the boundary, i.e., Ab = and Qb = Qo, we get back the equation (j5J). 

The above expression (j!2p simplifies dramatically if the operator O is located 
exactly on the boundary and if, at the same time, one can find a Laurant matrix A b 
satisfying 

Q'-Q b = A' b 77, 

with a projection rj that maps into some subspace of j = ojjj The equation (ITUj) then 
becomes simply 



Z[Q>] 



n 



sinh 2K 



j,k 



Z W \ it sinh 2K b jtk 



det 



n . 

i + Y- K V 



(13) 



A similar result is obtained for Q by simply putting on every matrix and reversing the 
order of matrix multiplication. 



2.3. Boundary operators in CFT 

We briefly summarize the CFT results on the boundary correlation functions of bcc 
operators appearing in the Ising model. At the critical temperature the continuum 
limit of the Ising model is described by the rational CFT of central charge |. There are 
three primary fields: the identity, the energy, and the spin operator with (holomorphic) 
dimensions 0, ~ and -h respectively. And there are three conformally invariant boundary 
conditions of which the two are identified with the two fixed boundary states, and the 
other is identified with the free boundary state. 

We summarize the dimensions of bcc operators and their correlation functions |16j . 

fixed-to-fixed : h+- = ~ (00) ~ \x\ — , 

fixed-to-free : h±f = i (00) ~ \x\ — x 2 \~* , (14) 

where X\ and x 2 are the points on the real axis (boundary) where the bcc operators are 
located. 

Let us remark on a subtlety when performing a lattice calculation of these 
correlation functions. Usually a correlation function (00) is determined by the ratio 
Z' jZ where Z' differs from Z by an insertion of 00. When calculating a correlation 
function such as the spin-spin correlation function we separate the effect of 00 from 
that of the background (the bulk and the boundary) by taking the ratio Z' jZ and thus 



§ This is possible for most of boundary correlators but may not be so, for instance, for four point 
functions of mixed operators. 
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Figure 1. Two lattice realizations of (fixed | free | fixed). The physical side of the half- 
plane is shadowed (in yellow). The thick bond is for infinitely strong coupling and the 
dashed bond is for vanishing coupling. The numbers below indicate j coordinates. A 
free boundary is inserted (a) by altering the vertical bonds and (b) by altering the 
horizontal bonds. Both realizations result in the same configurations in the physical 
(shadowed) region and, therefore, the same correlation functions. The length of the 
inserted boundary (free boundary) is n — k 2 — k\ + 1 which, in the case illustrated, 
equals four. 



cancelling out the bulk and the boundary free energy. Without proper cancellation the 
correlation function will contain a contribution that behaves exponentially with the size 
of the system. 

For the boundary correlation function of bcc operators (<f)<f)), however, the full 
cancellation of the background may not occur since the added operators can change the 
boundary state and thereby alter the boundary free energy. In this case the simple ratio 
Z'/Z will have an exponential behavior coming from the difference of the boundary free 
energies in Z' and Z. To extract the correct correlation function we must divide away 
the exponential factor that comes from the boundary free energy difference. 

3. Boundary correlation function of fixed-to-free bcc operators 

Let us proceed to calculate the boundary correlation function of fixed-to-free bcc 
operators. By placing two such operators there are two possible situations: a finite free 
section inside a fixed boundary, and a finite fixed section inside a free boundary We 
will use the notations (fixed | free | fixed) and (free | fixed | free) to represent the expectation 
values for each case. We start with the former. 

Our starting point is equation ( {TBI) . From now on we will simply write p instead of 
p y , suppressing the subscript. 

We assign our lattice uniform but possibly anisotropic couplings: K for vertical 
and K for horizontal bonds. We also assign another coupling Kf to the free boundary. 
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3.1. (fixed\free\fixed) 

There are two ways to realize (fixed | free | fixed) on the lattice: one is to insert the free 
boundary by changing vertical couplings as in figure 1(a) and the other is by changing 
horizontal couplings as in figure 1(b). We first present the former. 

We already defined Qo for uniform couplings in ([6]). The fixed boundary may be 
imposed by setting 

K 0jk = e and K 0>k = -lne, (15) 

and by taking the limit e, e — > 0. Let us define Qb as the corresponding Q- matrix. 
A finite free section is inserted by changing the vertical couplings as 

K , k = Kf for h<k<k 2 , (16) 

and Q' is defined accordingly by implementing ( fTBT) into Qb- We set k 2 — k\ + 1 = n so 
that the number of the inserted free bonds is n. 

We define the correlation function using (|T3|) as 

Z[Q>] ^ 



(fixed | free | fixed) „ = lim* 



lim* 



fee**} 



lim* * (ee K i)- 2n det 



n 

Ab 



A' hV 



'ee K t)- 2n D n 



Ab X nA b 



lim* stands for the limit e, e — > 0. D n means the truncated Toeplitz determinant of size 
n (see Appendix). In the last equality we used the fact that the rank of the projection 
T] is In which comes from the number of the modified couplings when inserting the free 
boundary (2 comes from the r-space). 

The additional factor of (ee^)"™ comes from the following reason. The modified 
n vertical bonds (flEj) have the couplings K { in Z[Q'] and (—hie) in Z[Qb]; those bonds 
are weighted differently in the two partition functions. The additional factor is added 
to cancel this effect and compare the two configurations (one with the inserted free 
boundary and the other without) within the same set of couplings Kf. 

To evaluate the determinant we only need to know the leading terms of Ab and A' b 
in small (e, e)-expansions. They are given by 

1 (cosp — l)[(s + s cosp — c) s — c] \ 



A,, 



A' 



4ee 2 
1 



cosp + 



4ee 



?2 



'o~\ ship + cosp — 1)(1 — a 2) 



1 



H [cosh2Kf(l + 02 cosp — o"i sinp) — sinh2i^"f(o"2 + cosp — ia^ ship)]. 

2 £ 



The squared correlation function becomes 
(fixed I free | fixed) 2 n 



lim* (ee Kl ) 2 " D 



l+a 2 



cri- 



es cosp — cs 



s{c — s) snip 



+ 2i z U 



_ 2 cosh2iCf(l + cr 2 cosp — o\ sinp) — sinh2i^f(o" 2 + cosp — ia 3 sinp) 



1 — cosp + (£ 2 — 1) 1 / 2 (cosp — l)[(s + s cosp — c) s 
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cs cosp — cs + v^ 2 — 1 \ 1 



a 2 



s (c — s) smp 



ai- 



cosp + 1 



det 



e 1 (a 11 b 22 - auhi) 



(17) 



At the last equality we applied the projection (1 — o~ 2 )/2 to the term of the order 0(e 2 ) 
since its (1 + cr2)/2-projected complement will not contribute to the determinant in our 
limit. This procedure may be illustrated using the following simplified equation. 

e _1 \ / a u a 12 \ ( fin ftu 
J \ a 21 a 22 J V & 2i ^22 

where the crossed-out terms, b\\ and &12, do not contribute in the leading order, and 
the "crossing out" is equivalent to applying a projection ( ° ° J to the 6jj-matrix. This 
procedure will be useful again. 

After a simple unitary transformation we obtain the correlation function as 

(fixed | free | fixed) = J D n (a*) , 



where the symbol (see Appendix) a* is given by 



cs cos p — c s-\- 



— e 



-4K { 



cot 



s (c—s) smp 
1 



Let us present the other way of calculating the same correlation function: by 
changing horizontal couplings to insert free boundary as illustrated in[T](b). Here it 
is convenient to use the Q notation defined in (jlj). 

First we define Q by assuming uniform couplings. We add fixed boundary at j = 
by assigning 

Ko,k = K , k = -In I, (19) 

and taking the limit e — > 0. We define Qb as the corresponding Q-matrix. 

A finite free section is inserted by cutting off the horizontal bonds as in figure [H(b) 



K 



0,k 



e for ki < k < k 2 



(20) 



and taking the limit of e — > 0. Remember that we just changed n — 1 horizontal couplings 
(we set k 2 — k\ = n — 1) to insert n free bonds. We define Q' by adding these changes 
into Qb. In contrast to the previous setup the above procedure restricts the coupling of 
the free boundary to be Kf = K. 

Using (fT3|) the ratio of the partition functions is written ajjj] 



Z[Q>] 



— sinh 2e 



n-l 



det 



a, n 

l+r/A' b — 



Z[Q h ) ^sinh(21ne); 

where the projection r] is into k\ < k < k 2 and has the rank 2(n — 1). 

The modified n — 1 horizontal bonds (120]) have the infinitely strong couplings (— me) 
in Z[Qb] and vanishingly weak couplings e in Z[Q']; those bonds are not weighted equally 



The tilde " will be suppressed for the rest of this section though it is implicit in every matrix 
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in Z[Qb] and in Z[Q']. We multiply by an additional factor to remedy this effect and 
define the correlation function as 

Z[Q>) " 



(fixed | free | fixed) „ = lim* 



lim* 



fee 6 ) 1 " 1 x 



Z[Qi 



(21) 



(4e)«- 1 J D n _ 1 



. . n 



lim* stands for the limit e, e — > 0. The determinant is over the truncated matrix of size 
n — 1 which comes from the rank of the projection r). 

To evaluate the determinant we only need to know the leading terms of Ab and A' b 
in small (e, e)-expansions. 

i i 



A b 



s c — sc cosp \ 

ve-i J' 

CT 2 1 + 0~2 COSp — 0~\ Ship 1 1 + (72 COSp — 0~\ SYO.p 



2ee 2 2 2 e 2 2 

The leading term of A' b contains a projection (1 —02)j2 and we included the subleading 
term to evaluate the limit (otherwise the determinant will be zero). Performing 
the limiting process by substituting with these leading terms the squared correlation 
function (12TT) becomes 

(fixed | free | fixed) 2 n 



= lim*(4e) n - 1 J D n _ 1 [e" 1 ±fz kha gggjpi *j*2 IL + A _ I 



= ^-1 
= A 



1— <T2 2(l+(72 cosp— (7i sinp)n 
2 



L 1+22 

l+(5 2 -l)- 1 /2( S c_sc C osp) "r" 2 



CT2 cosp— Q-j sinp II 
2 At" 

2(1+(T2 cosp— (T! sinp)n 



-i , 2(1+(T2 cosp— (Ji sinp)II 

i + a 2 - cr 2 i +K 2_i)-i/2 (s c_ IccoS j,) 



l+(^ 2 — 1) 1 ' 2 (s c— sc cosp) 



In the second equality we multiplied the projection (1 + er 2 )/2 to the subleading term as 
its (1 — U2)/2-projected complement will not contribute to the determinant in the limit 
(see the explanation near the equation f|T7|) ). 
We obtain the correlation function as 



(fixed | free | fixed) n = y D n _i(a) 
where the symbol a is given by 

a(e ip ) = 1 + cr 2 - 0-2 



2(1 + cr 2 cosp — 0\ sinp) II 



1 + (£ 2 — 1) 1 / 2 (sc — sc cosp) 



(22) 



3.2. (free\fixed\free) 



Next we calculate the reversed boundary state (free | fixed | free), where a finite fixed 
section is inserted into the otherwise free boundary. Since the procedure is very similar 
to the previous section we will only describe the setup and the results. 

In figure [2](a) we insert the fixed boundary by changing horizontal couplings. We 
define Qb to incorporate the free boundary as 



K 



0,k 



e and K { 



0,k 



— lne 
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(a) (b) 

Figure 2. Two lattice realizations of (free|fixed|free). The physical side of the half- 
plane is shadowed (in yellow). The thick bond is for infinitely strong coupling and 
the dashed is bond for vanishing coupling. The numbers below are j coordinates. A 
fixed boundary is inserted (a) by altering the horizontal bonds and (b) by altering 
the vertical bonds. Both realizations result in the same configurations in the physical 
(shadowed) region and, therefore, the same correlation functions. The length of the 
inserted boundary (free boundary) is n = &2 — k\ + 1 which, in the case illustrated 
above, equals four. 



Q' is defined to incorporate the inserted fixed boundary by changing horizontal couplings 



as 



K 



0,k 



— In e for k\ < k < k 2 + 1 . 



Note that the number of modified bonds is n + 1 and the length of the inserted fixed 
boundary is n = k 2 — k\ + 1. Using (TT3"j) the correlation function is defined with the 
appropriate normalization as 



(free | fixed | free) n = lim" 



\ 



ee c 



\n+l 



z[Q'\ 
Z[Q h \ 



lim* 



\ 



(4e) 



-n-1 



Da 



+ 1 



After similar steps as in section 3.1 we get 

(free | fixed | free) „ = ^D n+1 {b), 
where the symbol b is given by 



6(e 



1 



cs cos p— s c+ \J £, 2 — 1 
s (c—s) sin p 



-COtf 



(23) 



In figure Wfo) we calculate the same quantity by changing the vertical bonds to 
insert the fixed boundary. We define Qb to incorporate the free boundary as 



K 



0,k 



e and K { 



O.k 



We define Q' to incorporate the inserted fixed boundary as 
Kq k = — In e for k\ < k < k 2 ■ 
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With an appropriate normalization factor the correlation function is obtained as 



(free | fixed | free ) n = lim* 



\T 1 ^[Qbj 

After similar steps to section 3.1 we get 



r^) n S = lim\ (ee K *) n D r 



(free | fixed | free) n = ^JD n (b*), 
where the symbol b* (e ip ) is given in a rather complicated form by 

b*(e-)= 1 -^ + e ~ Kt{ 't" ) x (24) 
v ; 2 2s sin 2 p v ' 

- y/£ 2 - 1 (cosh 2K { - sinh 2K { cosp)+cosh 2K { (eg — secosp)+sinh 2Kf (— secos p+sc— ss sin 2 p) 
5c cos p— sc+cosh 4/S"f (— sc cos p— ss)+sinh 4K{ (sc+ss cos p) 

x (— cr 2 cs — (Jxss sinp + sc)(l + o~i cosp + <7i sinp)(cosh Kf — cr 2 sinh Kf). 
4. Some identities and the Kramers- Wannier duality 



The symbols: a* ([18]), a ([22]), 6 ([23]) and 6* ([24]), that we have obtained in the previous 
section all depend on K and K, and a* and b* have further dependence on Kf, the 
coupling of the free boundary. By setting Kf = K we expect a* (and 6*) to produce the 
same correlation function as a (and b) since the latter has the boundary coupling of K. 
More precisely, we expect 



D, 



D n _!(a), D„(6*) =£)„+! (6). (25) 

Af =A 



Though the second identity is trivial from the physical point of view the first identity 
is not totally obvious. Its RHS corresponds to figure [Qb) and, there, the two separated 
fixed boundaries must be of the same sort since they are actually connected underneath 
the free boundary through the unphysical side. However, in LHS which corresponds to 
figure [T](a) the two separated fixed boundaries can be polarized in different spin states, 
i.e., spin-up for one side and spin-down for the other. The identity then means that the 
"up-down" configurations are of measure zero compared to the "up-up" or "down-down" 
configurations. 

This can be understood by considering a one point function of fixed-to-fixed bcc 
operator and the Kramers-Wannier duality. We will show that this one point function 
vanishes, which will mean that the "up-down" configuration is indeed of measure zero. 
The Kramers-Wannier duality says that the lattice with couplings {K} has the same 
partition function as the dual lattice with couplings {K*} = {arctanhe~ 2i<: }. If one 
changes K into K+i^ then its dual coupling changes as K* — > —K*, which is equivalent 
to exchanging the spin operators and the disorder operators [15]. Let us polarize the 
spins on the vertical boundary by imposing a positive infinite couplings K = oo along 
the boundary. Then we change the sign of one bond to K = — oo, which is to put a fixed- 
to-fixed bcc operator. In the dual lattice, the (vertical) bond with the infinitely strong 
coupling becomes disconnected (horizontal) bonds (i.e. K* = 0) and the (vertical) bond 
with the coupling K = — oo becomes the (horizontal) bond with the coupling K* = i|, 
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which is simply to put a pair of spin operators on that bond. The partition function 
then vanishes naturally because one point function of the spin operator vanishes. 

We can also understand the similarity of a* and b from the Kramers- Wannier 
duality. The similarity may be explicitly expressed as 

D n (a*) _ _ = 2 2n D n (b). (26) 

Kf =0, c<->c, s<-^s 

By imposing K{ = we have chosen that the inserted free boundary (four bonds in the 
figure 1(a)) have vanishing couplings, that is, disconnected bonds. Since the Kramers- 
Wannier duality exchanges the fixed bond with the disconnected bond one obtains the 
figure |2(a) as the dual configuration. For the other bonds, the duality acts asc^c and 
s <-> s. 

The above paragraph explain the equation (|26|) except the factor of 2 2n . To explain 
the factor let us explicitly write the Kramers- Wannier duality as follows. 

TJ e K «'i = Yl C osh K ( tanh K) lcngth (FK loops) (27) 

{s} (ij") bonds loops 

= 2 II e ^E e _2 ^ lcnsth * (Ising domain walls). (28) 

bonds* loops* 

For each case the partition function has two parts: the product over all the bonds 
and the summation over loop ensemble. Since our Toeplitz determinants are the ratio 
of two partition functions we may express them using the above form. If we express 
D n{ a *) \k{=o,— loops and D n (b) by the Ising domain walls then the contributions 

from the loop-summations are equal; only the contributions from the bond-products 
may differ. By noting that cosh pa |e Xo ° for — > oo and also that we normalized 
the correlation functions by dividing out the infinite factor e Koc , we can confirm the 
existence of the factor 2 2n . 

From the symbol a* we may extract more information than just the correlation 



function. If we series-expand J D n (a*) in powers of e~ iKt then each term collects the 
configurations with the same number of "boundary islands" which we will define as the 
region of oppositely oriented spins within the uniformly fixed boundary. The number of 
boundary islands is given by the power of e~ iKi . 

The maximal number of boundary islands corresponds to 1+( -~ 1 - > e~ 2nK{ JT n (a\ 



2 c y ± n\^l,2) 

where a* l2 is the (1, 2)-component of the symbol a*. This corresponds to the insertion 
of an alternating boundary within the fixed boundary. Therefore we get the following 
as a byproduct. 



(fixed|alternating|fixed)„ = 1+( 2 1} " yXfe)- (29) 

The boundary condition must alternate even times for the correlation function to be 
non-zero. 



5. Asymptotic behavior of the determinant 



The asymptotic behavior of a general block Toeplitz determinant is obtained by the 
generalized Szego's theorem which is proven by Widom [TJ] (see Appendix). This 
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theorem may not apply, however, when the (matrix) symbol contains some singularities. 
For a scalar Toeplitz matrix there is the Fisher-Hartwig theorem (see Appendix) which 
tells you the asymptotic behavior of the determinant when its symbol contains the 
Fisher-Hartwig singularities. For block Toeplitz matrices there seems to be no analogous 
theorem, which we will need to study the correlation function at the critical temperature. 

Interestingly, however, D(a*) can be transformed into a scalar Toeplitz determinant. 
Let us recall the definition of a* f|T8|) . 

*/ m ( 1 w(p)\ cs cosp-cs + y/C 2 - 1 

a ( e ) = „- 4 fc„„,P i > where w{p)- 



-e " 4Xf cot | 1 /' ' ' s (c — s) sinp 

By performing a simple similarity transformation we get 

D (a*) - det ( 4 Tn{w) )-det( 4 -^ Kl Uw)T n (cot f ) 

" ( j " l-e- 4 ^T n (cotl) I n )- det [l n I n 



det 



I n + e- 4Ki T n (w)T n (cot 



A further simplification can be made by using the fact: [T n (cot §)] = T n (tan \ ) which 
holds for an even n. Multiplying this matrix inside the determinant we get 

D n (a*) = (-l)*/ 2 D n (tanf + e" 4 ^w(p)), for even n. (30) 

We have used _D n (cot |) = (— l) n / 2 for an even n. 

Let us find out the asymptotic behavior of the scalar Toeplitz determinant (!30l of 
which the symbol is the scalar function, 

h fjP) = tan 2 + ^ csoosp-cs+y/e^l 
y ' 2 s (c—s) sinp 

According to the Szego's strong limit theorem (see Appendix) the asymptotic behavior 
of the determinant is given by 

D n (h) ~ G{h) n E(h) as n -> oo. 

The exponential factor is simply given by 



G(/i) = exp(log /i) = exp f ^ log h(e ip ) 

Jo Art 



10 

Let us evaluate (log/i)o assuming K = K and K{ = K for simplicity. 

(log/i) = / 27r ^log[2e _2X (c-cs + s-s cosp + dl . (31) 

JO 27T L V /J 

Near the critical temperature, X ^ K c — | ln(l + v2), this becomes 

(log fc) = |C G + log 2(V2 - 1) + J - K c ) log |K - K c \ + ■ ■ ■ , (32) 

where Cc ~ 0.915966 is called the Catalan's constant. 

This exponential factor, G(h) = exp(log/t)o, has the following clear physical 
meaning. Let us prepare the boundary condition such that all the spins to one side 
of a certain boundary bond are polarized into the up-state. Then the probability that 
the free spin nearest to the polarized spins will also polarize to up-state, is given by 
G(h)-\ 
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At the critical temperature, as explained near the end of the section 2.3, this 
exponential factor indicates the boundary free energy difference. In our case, (log/t)o is 
twice the boundary free energy difference (per bond) between the free boundary and the 
fixed boundary because our correlation functions are the square root of the determinant. 

Let us calculate the asymptotic behavior of the determinant D n (a*) ( 1301 at the 
critical temperature, K = K c = |arcsinhl. We also set the boundary couplings as 
= K c . Then we need to consider the following symbol of the scalar Toeplitz matrix. 

h c {e ip ) = a^=i sec f (yi - cosp + ^3 - cosp) . (33) 

One notices that this symbol has a jump discontinuity and can be decomposed as 

h c (t) = (-ty/ 2 b c (t), (34) 

where b c (t) does not have any Fisher-Hartwig type singularity (see Appendix) but 
contains a simple pole divergence at t — — 1. This pole singularity prohibits us from a 
simple Wiener-Hopf decomposition of the function \ogb c (t) and from using the following 
formula presented in Appendix. 

3 1 00 

D n (h c ) ~ e" (log/lc)o ^^exp [£ fc(log6 c ) Jfc (log6 c )_ fe ]6 + (l) 1 / 2 6_(l) 1 / 2 n a2 - /32 . (35) 

k=\ 

Instead, as we did naively in Appendix, let us try to evaluate the asymptotic behavior 
by using the Szego's theorem in a formal way. 

Let us consider the ratio D n (h c ) / D n (\/^i) . Its exponential factor, e^ oghc '°, was 
already obtained to be -Cq at (1321) . The constant factor is written as follows by formally 
applying the Szego's theorem. 

_ oo oo 

(f c °) n zgr^y = gx p [ £ k ^ h ^ l °z + £ i] • ( 36 ) 

fc=i fc=i 

We calculate (log h c )k for k ^ using integration by part as 
(log h c ) k = ^e~ ikp ( log Vi - log(l + e ip ) + log(Vl - cosp + ^3 ~ cosp) 

^ , ( _ i a l fc l 1. • H\— c2tt dp e~ ik P cos(p/2) 

= ~2k + ~~ k~ m + W 6 1Qt |fc|Sgn(A;) = J Q ^ V6 _ 2 X • 

In the above equation we considered log(l + e ip ) as lim e ^o+ log(l + e + e ip ) hoping that 
the final result may not depend on the specific limiting process. Substituting this into 
the equation fl36l) we can cancel out the diverging harmonic series and obtain a finite 
quantity as follows. 

oo _,2 oo / i\fc„ oo / i \fc 

(2 C \~ n Dn(h c ) _ a k v (-1) a k v (-1) i 

These sums can be evaluated using the relation J^^ p = i Y,kL-oo a \k\sgn-(k)e lkp and 
its integrated form, arctanh^|==- = J2 k ^o jkj elkp + as 



oo ^,2 



a t f dp f dp' . p — p' cos(p/2) cos(p'/2) 



fc=i 



^ jfe 7 27TJ 2vr e 



2 sin 



\/6 — 2 cosp ^6 — 2 cosp' 



0.03921085231, 



| ; ( Z l)^ = l arctanh l _C_o ~(-l£ = J og2 
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Figure 3. D n {h c ) / [(£ C G ) n 0.5506047747 x 71" 1 / 4 ] vs. n. It shows the convergence 
to one as n goes to oo. 



Collecting all these we finally obtain the asymptotic behavior of the determinant as 
D n (h c ) ~ (f C G ) n exp [0.03921085231 + ±arctanh^ - ^ - l -^-]D n (V^i) 

' 2 : C G ) n exp [0.03921085231 + ±arctanh^ - - fi^e^A- 3 )^ 1 ^ 
0.5506047747 x n" 1/4 , 



la 



G 



where A = exp[^ - C'(-l)] = 1.282427129... is the Glaisher-Kinkelin constant and C' is 
the derivative of the Riemann zeta function. Figure [3] confirms this asymptotic behavior 
of the determinant. 



6. Conclusion 



Starting from the formulation of |4j we obtained the general boundary correlation 
function in a 2 x 2 block Toeplitz determinant (TIB"]) . We then calculated the 
boundary correlation function of fixed-to-free bcc operator. Four lattice realizations 
of this correlation function lead to four block Toeplitz determinants. We explained 
some interesting properties between them using the Kramers- Wannier duality. These 
determinants could be transformed into a scalar determinant when the size of the matrix 
is even. Using this scalar Toeplitz determinant we calculated the asymptotic behavior of 
the correlation function at large distance, which agrees with the conformal field theory. 
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Appendix. Toeplitz matrix: Szego's theorem, and Fisher-Hartwig 
determinants 

Here we collect definitions and theorems about Toeplitz matrices |17j . 
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Given a sequence {a n }^ = „ 00 of complex numbers let us define an infinite matrix M 
such that its elements are given by My = a^j for — oo < i,j < oo. Such matrices, 
that is, doubly-infinite matrices which are constant along the diagonals, are called 
Laurent matrices. Let us assume that there is a function a(t) defined on a unit circle, 
{t\t = e ip ,p G [0,2tt)} such that 



1 r2n 

— a(e ip ) e~ mp dp. 

2,71 Jo 



(37) 



The function a is referred to as the symbol of the matrix M. 

A Toeplitz matrix is an infinite sub-matrix of a Laurant matrix such as follows. 



T{a) 



clq a_\ 

ax a 

a 2 ai 

a 3 a 2 



a_ 2 a_ 3 

a_i a_2 

CLq Q,—i 

CL\ CLq 



V 



\ 



(38) 



where the function a is defined by (1371) and is called the symbol of the matrix T(a) in 
this context. If the symbol a is a finite dimensional matrix function then the resulting 
matrix is a block Toeplitz matrix. 

We are mostly concerned with the truncated Toeplitz matrix of size n x n which will 
be denoted by T n (a). The determinant of T n (a) will be denoted by D n (a). 



det T n (a) = det 



«o 
a i 



a n-l 0>n-2 



a 



When the symbol a is sufficiently smooth, non-zero, and satisfies J 2n loga(e i; 
(i.e., zero winding number), the Szegd's strong limit theorem states that the determi- 
nant D n (a) has the asymptotic behavior given by 



D n (a) ~ G(a) n E(a) 



n 



oo, 



where G(a) and E(a) are defined by 

oo 

G{a) = exp(loga) , E(a) = exp A;(loga) fc (loga)_ A; , 



k=l 



using the fourier components, (loga) n := J 27T & loga(e ip ). This theorem is generalized 
to block Toeplitz matrices by Widom [14] . 

The symbols with pure Fisher- Hartwig singularities are characterized by two 
numbers a and (3 as 

w a Ae ip ) = |2-2cos(p-po)| a e i/3(p - 7r) , 

where a deals with the zero-modulus singularity and (3 deals with the jump singularity. 
For some range of a and /3, the asymptotic behavior of D(w a ^) is known to be 



DnM WT^) n 



n — > oo. 



(39) 
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G is the so-called Barnes G-function which is an entire function defined by 

oo 

G(z + 1) = ( 27r ) z /V^ +1 )/ 2 - c ^ 2 / 2 J] {(1 + ^) n e-^ 2/(2n) }, (40) 

n=l 

where C 7 ~ 0.57721... is the Euler-Mascheroni constant. 

In general, when the symbol contains certain singularities, the sum 
J2kLi k{loga)k(}oga)-k appearing in the definition of E(a) may diverge and the determi- 
nant may pick up a power-law behavior. Such symbol may decompose into the product 
of a finite number of pure Fisher-Hartwig singularities and a "nice" function (i.e., a 
function subject to the condition that applies to Szego's theorem). For an example, let 
us write 

a(t) = exp(log a) w p (t) b(t) w p (t) = {-tf , (41) 

and assume that b(t) is a "nice" function that does not have any singularity. Let us 
calculate the ratio, p n ff y by formally applying the Szego's theorem respectively on 
both quantities as follows. 

D (a) 00 00 

—H — - = G( a ) n exp [53fc(loga)jfc(loga)_ fc - ^ /c(log w /3 ) fc (logM7 /3 )_ fe 
^W) L fc =i k=i J 

oo oo oo 

= G(a) n exp ^(log&) fc (log6)_ fc + ^ k(logb) k (\ogwp)_ k + ^ k(logu7/j)fc(log &)_*], 
k=\ fc=i fc=i 

where all three terms converge since b is "nice." By using (logw^)^ = — f (1 — 8ko), we 
can simplify the above as 



oo oo oo 

= G(a)"exp [J2 A;(log6) fc (log6)_ fc + /?£( lo S% - /3E( lo S 6 )- 
fe=i fc=i fc=i 

oo 

G(a)" exp [ fc(log 6) fc (log 6)_ fe ] 6 + (l) /3 6_(1)^, 



fc=i 



where 6± comes from factorization, 6 = &+&_, where b + (t) = exp[Y^ < ^L (}ogb) n t n } and 
b-(t) = exp[J]^=_ 00 (log6) n t n ]. Combined with ( 1391) the above equation tells you the 
asymptotic behavior of the determinant, D n (a). This is only a physicist's way of 
presenting the more general theorem which is proved by Ehrhardt and Silbermann [19|. 
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